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I . INTRODUCTION 



A prediction interval is an interval that contains a 
future observation with a pre-specif ied probability. The 
limits of the interval are functions of known observations 
from a family of known distributions with given parameters . 
Though prediction intervals resemble confidence intervals , 
they differ from the latter conceptually. A confidence 
interval covers the value of a parameter of a distribution. 

A prediction interval, on the other hand, encloses the value 
of a random variable. 

A prediction interval for a single future observation 
and simultaneous intervals for independent and identically 
distributed samples from N(u,c 2 ) were derived by Hahn [2]. 

For correlated samples with a prescribed correlation struc- 
ture, prediction intervals were developed by Choi [3]. It 
was shown by him that many of Hahn's prediction intervals 
are valid even for the case where the samples were correlated 
with the specified correlation structure. 

In this thesis, the observations are assumed to belong to 
a multivariate normal family with mean vector 



p 

nxl 



(p , v , • • • ,y) T 



and correlation matrix 




2 



n-1 



( 1 * 1 ) 
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Let X (t= . . . -2 , -1 , 0 , 1 , 2 . . . ) be a stationary Gaussian 

markov process of order one with E(X t ) = y for all t and 

2 k 

the covariance between X and X = a p for all t and k. 
Stochastic processes of this type are used to model certain 
types of oceanographic and metereological phenomena and 
some aspects of the stock market. Any finite sub process 
will then belong to the above multivariate normal family. 

Chapter 2 deals with prediction intervals from a normal 
distribution and contains a summary of known results and 
a review of general linear hypothesis. In chapter 3 a 
prediction interval for a future observation in a first order 
markov process is derived and the properties of this interval 
are examined. The effect of various parameters such as 
variance, correlation and sample size on the prediction 
interval are also discussed in this chapter. Chapter 4 
contains a discussion of a simulation to generate serially 
correlated random variables, which are used to empirically 
verify the theoretical conclusions drawn. The results of 
the simulation are presented in chapter 5. Extensions of 
the results to a linear trend model are discussed in chapter 6. 
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II. PREDICTION INTERVALS FOR NORMAL DISTRIBUTION 



A. SUMMARY OF KNOWN RESULTS 

Consider a random sample of size N=n^+n^ drawn from a 

2 

normal distribution N(y,c ). The samples n^ form the group 

of initial samples and n^ form the group, called future 

_ 2 

samples. Based on the mean X and variance S of the 
r n. n. 

1 - 1 2 

initial samples, prediction intervals for and were 

obtained by Hickman [6]. Hahn [2] derived prediction 
intervals for the mean and variance of the second sample 
and also simultaneous prediction intervals for the variance 
of each of k additional random samples of size n^ based on 
the information obtained from the initial sample. Hahn's 
[2] two sided 100 r % prediction interval to contain a 
single additional observation X is given by 



Y ± t [n-1 ; (1+r) / 2] (l+\) h • S 
_ 2 

where Y is the sample mean, S is the sample variance and 
t[n-l , (1+r) / 2 ] is obtained from the tables of t distribution 
with (n-1) degrees of freedom. A 100 r % simultaneous 
prediction interval to contain X^, X 2 , ..... X,^, k future 

observations is given by 



Y + R(k,n,r)S 



where 

R(k,n,r) = U(l+^) % 

The value of the function U can be obtained from the special 
tables given in reference 2 . 
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A prediction interval for a single future observation as 
well as simultaneous intervals for a specified number of 
future observations have been derived by Choi [3] for the 
case where the observations are correlated and belong to a 
multivariate normal distribution with mean 

T 

y_ — » 

and covariance matrix, 

V - %( H + _H T ) + ct( I - E ) 
nxn nxn nxn nxn nxn 




h^(i=l, 2 , . . . ,n) and a are positive constants. E is an nxn 
matrix all of whose elements are unity. 

B. DEFINITIONS, NOTATIONS AND REVIEW OF SIMPLE LINEAR 
REGRESSION ANALYSIS 

A brief review of the results and notation in a simple 
linear regression model are presented below. These results 
are used in chapter 3 in deriving prediction intervals for 
a first order markov process. By way of notation, a 
capital letter with an underbar represents a vector or a 
matrix. Consider the simple linear regression model 

Y = X B + e (2-1) 

nxl nx2 2x1 nxl 

2 

where £ v N(0,a I) 
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which m scalar notation can be written as 



Y i = a + bX i + e. 



i = 1,2, ...,n. 
The least squares estimator for B is obtained by 



that is 



Min Q = (Y - XB ) 1 (Y - XB) , 
B 

~ -IT 
B = S 1 X i Y , 



where 



S = XX , 



and 



B ^ N(B, a 2 S _1 ) 



(2-2) 



B is the column vector (a,b) 

It is well known that the random variable 

Y ~ N (XB , a 2 I) 



and 



.-1 



n-SXX 



SX? -EX. 
i i 



•EX i 



n 



where 



n - 2 
SXX = E (X. - X) 

i=l 1 



2 

Let a be the usual unbiased estimator of o , given by 



a 2 = (SYY - b 2 SXX) / (n - 2) , 



(2-3) 



n 



V?\ 2 



where SYY = E (Y. - Y) 
i=l 1 



9 



A 

Let Y be the predicted value of the random variable at a 
specified value of the independent variable x. 

A 100(l-a)% prediction interval for the dependent vari- 
able Y can be obtained from the relation 



P 




1 - a (2-4) 



Equation (2-4) will be used in deriving a prediction interval 
for x in chapter 3 . 

C. STATIONARY GAUSSIAN MARKOV PROCESSES 



Let X t (t = . . . -2 , -1 , 0 , 1 , 2 , . . . ) be a stationary gaussian 
markov process, i.e., an autoregressive stochastic process, 
of order one with E(X t ) = m and covariance between X t and 



X 



2 k 



t+k = a p • 

To derive the prediction interval for a single future 
observation ^k+l based on t ^ ie observations up to the 

results of Ogawara [1] will be used. He has shown that 
when the values of X 2 ^_^ (k=l , 2 , . . ,n+l) are fixed, X 2 ^ 

(k=l, 2, . . . ,n) are mutually independent. The conditional 
probability densities are given by 



f(X 



2k 



X 2k-1 ’ X 2k+1 ) 



V2 



EXP 



it a 



where 

a = m(l-p 2 ) / (1+p 2 ) 

b = 2p/(l+p 2 ) 





(a+bxp 



2 



(2-5) 
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a 2 - o 2 (l-p 2 )/ (1+p 2 ) 

X k ■ < X 2k-l + X 2k + l)/ 2 



A prediction interval for be derived in chapter 

3 based on this conditional distribution of the even numbered 
observations conditioned on the odd numbered ones. 
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III. PREDICTION INTERVAL FOR A FUTURE 
OBSERVATION IN SERIALLY CORRELATED NORMAL SAMPLES 



A prediction interval to contain a single future obser- 
vation for samples from a multivariate normal distribution 
with a correlation matrix of the form defined in equation 
(1-1), is derived in this chapter. In section A a conditional 
prediction interval to contain a single additional observation 
based on the conditional distribution of X 2k is obtained. 
Section B deals with the properties of the prediction interval, 
In section C the dependence of the prediction interval on 
sample size, p and sigma are discussed. 



A. DERIVATION OF THE PREDICTION INTERVAL 

Assume a sample of observations, X-^, X 2 , X^ , ..., ^k-l’ 

X 01 from a normal distribution with mean zero and covariance 
2k 

1c 

matrix V such that Cov (X t , X t+k ) = p . 

The conditional probability densities, as indicated in 
equation (2*5), are 



f(X 2kl X 2k-1 ’ X 2k+1 ) 



>/ 2ttc 



EXP 



yy? { x 2k‘ < a+bX P 



with conditional mean 



E(X 2k| X k ) = a+bX k ’ 



and X^ = (X 2k . 1 + X 2k _ 1 )/2 



(3-1) 



E (X 



2k 



X^) considered as a function of X k is callled the 



regression function of X 2k on X^. Graphically, it represents 
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the locus of the center of gravity of the conditional random 

*k- 

As a result of Ogawara's Theorem the following conditional 
simple linear regression model may be assumed, 

X 2k = a + bX£ + e k , k = 1 , 2 , . . . , n (3*2) 

2 

where s k = N(0, a ) . 

The maximum likelihood estimators of the parameters a, b 
and a 2 are given by 

A 

b = SXY/SXX 

a = x 2k " ^ X k 

a 2 = (SYY - b 2 SXX) / (n-2) , 




where 

— , _ 1 / x l~ t ~ x 3 + x 3 +x 5 + , x 2k-l +x 2k+l 

x k n \ 2 2 ’ ’ ’ ' 2 



x 2k n ( x 2 + x 4 + " • + X 2k ) 





n 


/ \2 




2 




sxx = 


^ i 

k=l 


( X 2k ) 


- n 


x k 






n 










SXY = 


E 

k=l 


X k X 2k - 


n x/ 
k 


x 2k ’ 




SYY = 


n 

Z 

k=l 


( X -) 2 


- n 


- 2 
x 2k 


(3-3) 



At a known value of x k the predictor x 2k is given by 

x 2k = a + b x k (3-4) 
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A 100(l-a)% prediction interval for the dependent variable 



x 2k can obtained from the relation 



x 2k “ x 2k 



- / i 

°V 1 + K + — 



- x k ) 2 



< t 



sxx 



= 1 - a . (3*5) 



From the above relation a prediction interval for x^.^ can 
be obtained noting that x^ is a linear function of * 

The inequality inside the brackets can be written as 



(x 2k ' x 2k ) 5 C 0 






+ 



( x k - x k) 2 
SXX 



(3-6) 



Substituting for x^ = (x^ ~ b x£) + b x£ and squaring both 
sides, the inequality (3*6) reduces to 



(x 2k ~ x 2k^ ^ ^ Xv ” x v) ” 2b (x 9V - x 9v ) (x^ - x^) 



k k 



2k 2k' v "k k' 



< t 2 a 2 (l + i) + t 2 a 2 (x£ - x^) 2 /SXX 



or 



SXX / W k 



- 2b < x 2k - x 2k )(x k - V 



(x 2k ' X 2k )2 



t 2 o 2 (l + 1) < 0 (3-7) 



This is a quadratic equation of the form 



acx: - x/) 2 + b(x’ - <) + c < o 



(3. 8) 
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where 



A = 




v 2~2 
t a 

SXX 



B = _2b (x 2k - x 2k> ' 



r _ - n 2 2~2 n+1 

C " (x 2k ' x 2k ) - c a — 



(3-9) 



Thus, the probability statement in equation (3-5) is equiva- 
lent to 

P[A(x^ - 5^) 2 + B(x£ - xp + C < 0] = 1 - a . (3*10) 

A prediction interval for x£ can now be obtained in terms of 
the roots of the quadratic equation (3 - 8) viz., 



K - 



-B±/ B 2 - 



4AC 



2A 



(3-11) 



the term under the radical sign simplifies to 



2 2~2 
B - 4AC = 4t a 



x 2k x 2k) + n+1 f £2 t 2 a 2 



SXX 



n 



SXX 




(3-12) 



Since = (x 2k . 1 + * 2 k+l )/2 , 

The unknown future value X 2 k+]_ can be written in terms of the 
known value as 



*2k+l = D ± E w ^ ere 
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(3-13) 



E = 2t c 



/ 




(3-14) 



Thus a prediction interval for x 



2k+l 



is given by 



D - E < x 2k+l - D + E 



(3-15) 



where D and E are given by equations (3 ’13) and (3-14) respec- 
tively. The properties of this interval are discussed in 
section B. 

B. PROPERTIES OF THE PREDICTION INTERVAL 

The prediction interval in the previous section was 
obtained by solvxng 



where the coefficients A, B and C are defined in (3-9). The 
roots of the equation and hence the form of the prediction 
interval will depend on the relative sizes of A, B and C. 

The effect of these coefficients on the prediction 
interval will now be examined breaking in to four all in- 
clusive cases. 

Case 1. A > 0 and B^ - 4AC > 0 



A(x k ~ K )2 + B(X k - + C = 0 




16 



or b 2 SXX 

2 A 2 
t a 



> 1 



^2 

T _ b sXX _ „ 

Let — Tps: — j — F 

t a 



(3-16) 



B > 4AC implies that 



4 b 2 (x 2k - x 2k )2 > 4 (V 



. 2 ~ 2 

t a _ 2 

SXX } (x 2k " x 2k ) 



4 (b‘ 



2 - 2 2 ~ 2 
t a t a 



SXX y n 



(n+1) 



2 ~ 2 \ ^ 2 ~ 2 



or \ b‘ 



SXX 



n 



(n+1) > - 



. 2 ^ 2 
t a 



( x oi, ~ X 0 , ) 



SXX v 2k 2k' 



~2 

b : s F - 1 > 


n 


CM 

>< 


- x 2k> 2 




n+1 


2 


T~ 


t a 


t 


o 


or i \5 XX > 


1 - n 


(x 2k 


' X 2k )2 


or 2" 2 


1 n+1 


~T 


A 2 


t a 


t 


0 


To satisfy the conditions A 


> 0 


and B 2 : 


that 


F > 1 and F > 1 


n 


CM 

X 


CM 

a: 

CM 

IX 

1 


n+1 




^2 

a 


Since the quantity 


(x 2k • x 


N 2 , . 2~2 . 
2k ) /to is 



(3 • 16a) 



conditions in equations (3- 16a) are equivalent to the condition 
F > 1. In this case the prediction interval for x 2k+ ^ will 
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be of the form 



[(D - E) , (D + E)] . 



This type of two-sided interval is referred to as TYPE 1, 
where D and E are as defined in (3 * 13 ) and (3’14) respectively. 



Case 2, 



a - 4AC > 0 and A < 0 



A < 0 



. 2 ~ 2 
t a 

SXX 



< 0 



or b < 



. 2 ~ 2 
t a 

SXX 



or 



^2 

b sXX 

T1 

t a 



1 . 



i . e . , A < 0 



F < 1 



(3-17) 



2 n (x 2k ' X 2k )2 

a - 4AC > 0 =*• F > 1 - -^r — — 

n+1 t Z a 2 



(3-18) 



Thus, if F < 1 and F i 1 ■ 



(x 01 - x 01 ) Z 

n 2k 2k 



n+1 TT 

t a 



the resulting interval will be of the form, 



(- «, D + E] [D - E, oo ) 



(3*19) 



and will be called a TYPE 2 interval. 



Case 3. If B Z - 4AC < 0 and A 0 the interval is called 



TYPE 3 
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The roots of the quadratic are complex and the prediction 
interval must be taken as (- ~, + «>) . This situation occurs 



when B - 4AC < 0 . 



B - 4AC <0 * F < 1 



n (x 2k' x 2k )2 



n+1 



. 2 A 2 
t o 



(3-20) 



Thus a type 3 interval results when 



F < 1 - 



n (x 2k ' x 2k )2 



n+1 



JT 2 
t a 



CASE 4. 



A = 0 
A = 0 



~2 

b SXX 

~7Tl 
t a 



= 1 . 



(3.21) 



Thus, when F = 1, the interval is of type 4 and is of the 
form [G,°° ) , 



where G = 2x ’ - x 9 , - b SXX (— ) + * 2k - A - X ? k 

k 2k- 1 n , 



(3-22) 



To summarize the results, the quadratic equation that 
leads to the required prediction interval is of the form 

A(x k " x k )2 + B(x k - x£) + C = 0, 



and four different types of intervals can result depending 
on the value of 



F = 



~2 

b z SXX 
t a 



The four cases are 
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i) 

ii) 

iii) 



iv) 



when F > 1, the interval is of the form 



C(D-E) , (D+E)]. 
when F < 1 and F > 1- 



n (x 2k ‘ x 2k )2 



n+1 



~T r T 

t a 



, the interval 



is of type 2 and is of the form (- EH-E] , [D-E , <») 



when F < 1- 



n (x 2k - x 2k> 2 



n+1 



. 2 A 2 

t a 



, the interval is of 



type 3 and is of the form (-<*>, +») 

when F = 1, the interval is of type 4 i.e., [G,«>) . 



At a given level of significance a, from the relation 
(3-5) and the conclusions drawn in this section, the following 
relation should hold: 

4 

Z Pr [occurrence of type i interval] 
i=l th 

x Pr [i interval contains the predicted value] 

= 1 - a (3-23) 



The above conclusions are verified with simulated 
samples in chapter 5. 

Section C examines the distribution of the quantity F 
defined in equation (3*16). 

C. DISTRIBUTION OF 'F' 



The occurrence of any particular type of interval depends 
on the value 



F = 



A 2 

b z SXX 
1~T~ 

t a 



The distribution of F can be shown to be a noncentral F 
distribution as follows. 
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From (2-2) , 

A 



or 



(b-b) 



or 



/v ry 

(b-b) ^ 
cr 2 /SXX 



. A 2 

it can be shown that b is N(b,a /SXX,) 



* N(0 , 1) 

v x 2 d), 



b 2 - SXX „ 2 ,, 

and j ^ X (1, A) 

a 



where 



b 2 SXX 

2 

a 



(3-24) 



is the noncentrality parameter. 

From (2-3) the distribution of 
~2 

2 ^ -— is X 2 (ji-2) . (3*25) 

a " 1 



The ratio of (3-24) and (3-25) devided by their respective 
degress of freedom results in a noncentral F [5] statistic 
with degrees of freedom (1, n-2) and noncentrality parameter 
X . That is 



~2 

b • SXX 




t 2 F ^ F (1, n-2, X) 



(3-26) 



and the noncentrality parameter 

, _ b 2 SXX 
A ^ • 

a 

Substituting for b from equation (2 * 5 ) , the noncentrality 
parameter X can be written as 
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X 



(3-27) 



4p‘ 



d+p^r 



sxx 

a 



D. NUMERICAL EXAMPLE 

The Model developed in section A is tested on the data 
obtained from Dow-Jones monthly averages to predict a future 
monthly average, based on the averages of the previous months. 
In most of the cases, the future value was contained in the 
prediction interval. To illustrate an example the data for 
the years 1966 and 67 are given below. The prediction 
intervals are computed with the help of the APL program 
PREDICT, whose listing is given on page 68. 



1966 1967 



Month Ended 


D - J Ind 


Month Ended 


D - J Ind 


January 


31 


983.51 


January 


31 


849.87 


F ebruary 


28 


951.89 


February 


28 


839.37 


March 


31 


924.77 


March 


31 


865.98 


April 


29 


933.68 


April 


28 


897.05 


May 


31 


884.07 


May 


31 


852.56 


June 


30 


870.10 


June 


30 


860.26 


July 


29 


847.38 


July 


31 


904.24 


Augus t 


31 


788.41 


Augus t 


31 


901.29 


Sept 


30 


774.22 


Sept 


29 


926.66 


October 


31 


807.07 


October 


31 


879.74 


November 


30 


791.59 


November 


30 


875.81 


December 


30 


785.69 


December 


29 


905.11 
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The above data are stored in the variable DMA. The pro- 
gram Predict is run with different sample sizes and the out- 
put is shown on pages 24 and 25. 

Using 14 monthly averages from January ' 66 to February 
'67, prediction intervals tor March '67 are computed and the 
exact value of 865.98 is contained in the interval. 

Using given values up to April ‘67, the prediction in- 
tervals for May '67, also contain the true value. With a 
sample of size 18, the prediction intervals for the nine- 
teenth month i.e., July '67 are computed and they contain the 
true value. Though the average of September '67 falls outside 
the interval, the average of November '67 is contained within 
the interval. 

The reduction in the length of prediction interval with 
increase in sample size, is supported by the above example. 

The data are tested for the required correlation structure 
and found to fulfill the requirements partially. The statistics 
of the data are 



Mean 



870.83 



Standard deviation 



55.0 



Corr. Coefficient p 



0.82 



23 



PREDICT 14tDMA 



THE F STATISTIC IS 
THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 



4.132073877 

727.1128824 

1080.631201 

353.5183188 

865.98 



PREDICT 16+DMA 



THE F STATISTIC IS 
THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 



5.489038103 

598.3935849 

900.5564036 

302.1628187 

852.56 



PREDICT 18 + DMA 



THE F STATISTIC IS 
THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 



5.130517543 

708.7039 

1009.101896 

300.3979963 

904.24 



PREDICT 20 + DMA 



THE F STATISTIC IS 
THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 



5.777973921 

573.5929826 

864.8705316 

291.277549 

926.66 
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PREDICT 22-t-DMA 



THE F STATISTIC IS 
THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 



6.994239226 

633.1176654 

899.7280383 

266.6103729 

875.81 
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IV. SIMULATION OF SERIALLY CORRELATED SAMPLES 



To study the effect of p , o and n on the type of pre- 
diction interval obtained it is necessary to generate random 
variables with the required variance covariance matrix. 

These can be generated as linear combinations of standard 
normal random variates . 

A. SIMULATION FROM STANDARD NORMAL SAMPLES 

Suppose Zp Z 2 , • • • , Z n are random samples from N(0, 1). 
To get serially correlated random variates X-j , , • • • >X n , 

let X 1 






X n = C i Z l +C 2 Z 2 + 

The coefficients C-^ , . ..,0^ 
equations of the form, 



+C!Z . (4-1) 

n n 

are then determined by solving 



Cbv (X t . X t+k ) 




for all t and k from 1 to n. 



Consider the first equation 
X 1 = Z 1 

E(X 1 ) = ECZ-l) = 0 

Var (X 1 )= Var (Z x ) =1 (4' 2) 
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From the second equation 



E(X 2 ) = E(C 1 Z 1 +C 2 Z 2 ) 

= c 1 E(Z 1 ;+c 2 e(z 2 ) = 0 

since Z^ and Z 2 are independent; 



Var (X 2 > = Var(C 1 Z 1 + C 2 Z 2 ) 

2 2 

— E ^ V ar Z x i" C 2 V ar Z 2 ( i . e . , 

2 2 
C 1 + c 2 = 1 



Cov(X ] _, X 2 ) 



E C(x 1 - E(X 1 ))(X 2 - E(X 2 ))] 
E [X 1 - X 2 ] 

E [Z 1 -(C 1 Z 1 +C 2 Z 2 >] 

E [E^Z^ +C 2 Z^Z 2 J ; i.e, 



C 1 E(Z 1 2 )+C 2 E(Z 1 Z 2 ) = p . 

Since * N(0, 1) , Z ± 2 v X 2 (l) ; 



(4-3) 



(4-4) 



also since Z^ and Z 2 are independent, the above equation can 
be simplified to 

C 1 = p (4-5) 

From 4.4 by substitution we get C 2 = yl - p 2 



or X 2 



pZ n + y 1 - p' 



(4-6) 



or X 2 




(4-7) 



27 



Similarly it can be shown that 



X 



3 




(4-8) 



The result can be generalised in the form 

x i = px i-i + V 1 ■ p2 z i < 4,9 ) 

i = 2, . . . , n. 

or X ± = p 1 " 1 Z l + p 1_2 y 1 - p 2 Z 2 + + {l - p 2 Z j _ (4-10) 

i = 1, . . . , n. 

The results obtained in (4-10) can be expressed in the follow- 
ing matrix notation : 

X = A Z , (4-11) 

nxl nxn nxl 

where X = (X, , X 0 X) T , 

nxl 1 z n 

Z — (Z, ’ ^9 * ' ’ ’ ’ ’ 

nxl 



and A is a nxn matrix of the form 




If a random variable Z is distributed normally with 

nxl 

mean jj and covariance V and if A is an nxn matrix of rank n, 
then the vector 
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Y « A 
nxl nxn 



Z 

nxl 



is distributed as multivariate normal with mean vector A y 
and covariance A V A . As a result, the random variable X 
(4-11) is distributed as multivariate normal with the re- 
quired covariance matrix V. Here 

Z v N (0 , I) , 

where I is an nxn identity matrix. 

Hence, X v N(0, A I A T ) . 

T T 

We now show that A I A = A A is a matrix of the required 
form (for n = 4) 



A 

4x4 



0 






p /i-p 

p2 p 



0 

0 



/l-p 2 /l-p 2 



p 2 /l-p 2 p/l-p 2 



0 

0 

0 

/l-p" 



A T - 

4x4 



1 

0 

0 

0 



/l-p 2 p/l-p 2 p 2 /l-p 2 
0 /l-p 2 p/l-p 2 

0 0 /.-p 2 
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1 p 



A A T = 



p p 2 +l-p 2 
p 2 p 3 +p(l-p 2 ) 
P 3 pVp 2 (l-p 2 ) 



p 3 +p (,1-p 2 ) 
p 4 +p 2 (l-p 2 )+(l-p 2 ) 
p 3 +p 3 (l-p 2 )+p(l-p 2 ) 



4, 2,, 2 S 

p +p (1-p ) 

p 5 +p 3 (l-p 2 )+p(l-p 2 ) 
p^+Cl-p 2 ) (p 4 +p 2 +l) 



which reduces to 



A A T = 



P 

1 

P 



P 

P 

1 

P 



P 

P 

1 



The above result can be generalized for any value of n, 



B. TESTING FOR NORMALITY AND CORRELATION 

The simulated samples are tested for normality and the 
correlation structure. The simulated samples are tested by 
standard Kolmogorov-Smirnov test. The results showed that 
samples do appear to be from the specified normal distribution. 
The correlation between and Xj is given by 

Cov(X. , X.) 

R = 1 ■ (4-13) 

VVar X i • Var X. 

Confidence limits for correlation coefficient p are developed 
using the approximation. 
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Z - % in i±| * H(% In jlj) 



A y57o confidence interval for p is given by 



U 2a -l)/(e 2a +l) < P < U 2b -l)/(e 2b +l), 



(4-14) 



where 



a = Z - and b = Z + 



/n-3 



/n-3 



The simulated samples are also tested for the required 
correlation structure by examining the confidence limits for 
correlation between (X^, X ^ ) using equation (4-14). The 
results showed that the samples do appear to have the pre- 
scribed correlation structure. 
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V. RESULTS OF SIMULATION 



In this chapter, using simulated samples, the probabili- 
ties of occurrence of various types of intervals and their 
dependence on the value of p , n and o are studied. The pro- 
gram PROBS FORTRAN, whose listing is given on pages 63-66 
computes the probabilities for the four types of intervals. 

For given values of p , n and a the program generates one 
thousand samples and records the frequency of occurrence of 
each type of interval. In each case the program tests whether 
the predicted value is contained within the interval. Finally, 
for a given level of significance a the relation (3 * 23) is 
verified . 

Choosing a level of significance a = 0.05, the simulation 
was run for values of p ranging from 0.1 to 0.9, for sample 
sizes from 10 to 50 and for sigma varying from 1 to 5 . For 
a = 0.1 and 0.02 the program is executed to study the effect 
of the level of significance on the results. 

The computer output is shown on pages 43 to 57 . Tables 
I to VI show the variation of probabilities of types of 
intervals with change in a at different sample sizes, keeping 
a and a constant. Tables VII to IX show the variation of 
probabilities of types of intervals with a change in p at 
values of a = 0.1. Tables X to XII show the variation of 
probabilities of types of intervals with a- change in p , 
keeping a = 0.02. 
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In Tables XIII to XV, the variation of the probabilities 
of types of intervals are shown at a standard deviation (a) 
of 3.0, whereas in Tables XVI to XVIII, the value of standard 
deviation is altered to 5.0. Tables XIX to XXVII depict the 
variation of the types of intervals with change in sample 
sizes at different values of p and o, keeping a at 0.05. 
Finally, Tables XXVIII to XXX show the variation of proba- 
bility of types of intervals with change in sigma at different 
values of p and n. 

In all these tables the first column indicates the value 
of the parameter. Columns 2, 3 and 4 give the probabilities 
of occurrence of types of intervals. Columns 5 and 6 show the 
probability that the predicted value is contained in the 
interval. Since a type 3 interval is of the form 
the probability of X 2 ^ + -i to be contained in the interval is 
one. Finally column 7 verifies the relation (3.23). The 
attached graphs on pages 58 to 62 show the trend in the pro- 
bability of occurence of the three types of intervals with 
changes in parameters like p, a, n and a. 

The outcome of the simulation confirmed the following 
theoretical conclusions : 

i) The probability of type 1 interval increases with an 
increase in correlation coefficient p, all other parameters 
remaining constant. 

ii) The probability of type 1 interval increases with an 
increase in sample size, keeping all other parameters constant. 
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iii) The probabilities of types 2 and 3 decrease with 
an increase in correlation coefficient or sample size. 

iv) The probabilities for the four types of intervals 
do not change significantly with change in sigma. 

v) The probability of a type 1 interval increases with 
an increase in level of significance a, all other parameters 
remaining constant. 

Section 6 deals with the application to linear trend 
model . 
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VI. APPLICATION TO LINEAR TREND MODEL 



In this chapter, the results of chapter 3 are extended 
to stochastic processes that follow a linear trend model. 
Let {X^} be a discrete stochastic process satisfying 



X 



n 



= y + e 
n n 



where y = a + Bn , 
n 



( 6 - 1 ) 

(6-2) 



and e n is a stationary process of normal variates that 
satisfy the following: 



E{e n> - 0 

Var{e } = c 2 
n 

Cor{e n> e n+k } = a 2 p k , p < 1 . (6-3) 

Such a process is called a model for linear trend. Krishnaih 
and Mur thy [7] showed that for the above model the conditional 
distributions of X 2k given X 2k-1 and X 2k+1 are independent 
normal with 



E [X 2k ! X 2k _ 1; X 2k+1 J B l + 3 2 K + B 3 X k 
and Var [X 2k | X 2k _ 1 , X 2k+1 ] = a 2 (1-p 2 ) / (1+p 2 ) , 
where X£ = (X 2k _ 1 + x 2k+1 >/2 



(6-4) 

(6-5) 
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In section A, a prediction interval for X 2k+ -^ 
obtained knowing the values up to X 2k for the above model 



A. PREDICTION INTERVAL FOR X 2k+1 

Using regression analysis, a prediction interval for 
X 2 k + i can be obtained for the linear trend model, knowing 
the values up to X 2k - Equation (6-4) can be written as 



X 2 k = 61 + 3 2 k + B 3 X- + E k , k = 1 , 2 , . . . , n, 



or in matrix notation, 



X ryt = X B + £ , £ ^ N(0 , a 2 I) . ( 6 - 6 ) 

nxi nx3 3x1 nxl 



The design matrix X is of the form: 

nx3 

1 1 (X 1 +X 3 )/2 

1 2 (X 3 +X 5 )/2 

1 3 (X 5 +X y )/2 



1 n (X 2 n . 1 + X 2n+1 >/2 



T 

and B is the column vector ( 8 ^, 8 2 , 6 3 ) 

3 x 1 

A 

The estimator B can be obtained from the normal equations 
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Pp A 

X X B = 
3xn nx3 3x1 



* T — 2k 
3xn nxl z 



or 



B = 



-1 T 
S A X A X 



2k 



where S = 



3xn nx3 





n 


Zn 


ZX’ 

n 


= 


Zn 


Zn 2 


ZnX’ 

n 




EX' 


ZnX ' 


EX’ 2 




n 


n 


n 



(6-7) 



Since S can be shown to be nonsingular and of rank 3, S - '*' 
exists and 

B v N (B, a 2 S -1 ) . (6-8) 



For a given value of X', 

1x3 K 



= [1 k x£ ] , 



the predicted value X = X', B 

lxl ZK 1x3 K 3x1 



and X 2k ^ N {X k B , a 2 (X^ )S _1 (X^ ) T } 



It can be shown that the residual 



x 2k " x 2k ^ N{0 ’ a2(1+ s 2) } > 



(6-9) 



where 



s 2 = X’ S _1 ( X' k ) T 
lxl lx3 X 3x3 3xl K 



As an example, let X = (x-^, x 2 , • • • » x ]_q) 

Here K = 5 and n = 4. 



The vector X 2k = ^ x 2 ’ x 4’ x 6’ x 8^ 

4x1 z 
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1 



and Che design matrix X = 

4x3 



1 

1 



1 



1 (x 1 +x 3 )/2 

2 (x 3 +x 5 )/2 

3 (x^+Xy) / 2 

4 (x 7 +x g )/2 



XjJ. =[15 (x g 4x 11 )/2] 



The value of x^ is unknown. A prediction interval for x^ 
can be developed as shown below. 



From equation (6-9) 



X 2k " x 2k 
a/ 1 + s 2 



N (0, 1) 



( 6 - 10 ) 



and 




t . 



( 6 - 11 ) 



The quantity t has a student's t distribution with (n-3) 

2 

degrees of freedom, where is the usual unbiased estimator 
of a 2 . 

At a given level of significance a , 





U 2k - x 2k 5 ■ /n - 3 


< t /0 
a/2 




o 




R / 1 + s Z 
o 








__ 



1 - a (6-12) 



The term in the brackets of the left hand side of equa- 
tion (6-12) can be written as, 
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(6-13) 



x 2k ' x 2k 4 Va/2 



and 



1 + s 2 = 1 + [1 k x^][S _1 ][l k x^] T . 



The expression on the right hand side can be simplified to 
the form 

Y o + Y l x k + Y 2 x k 2 ' 

The coefficients y , y^ and y^ can be computed since k and 
S ^ are known. 



Squaring both sides of equation (6-13) and substituting for 

A 2 

x 2k and s , it can be simplified to. 



t 2 R 2 

(x 2k‘ 6 l B 2 k_B 3 x k ) ^ n-3 ^ Y o +Y l x k +Y 2 x k ^ 



A A A 



or 



t - 2 r 2 t 2 r 2 y 

2 v a/2% „.2 . r2fx _g ~ k) : . V2V1, , 

(B 3 " y 2 n-3 } *k + (2(x 2k B 1 B 2 k;B 3 n-3 ' \ 



2 2 

tf/ofCY, 



+ < < x 2k -S 1 -6 2 k) 2 - > 



« 0 . 



The above inequality reduces to 



A x^ 2 + Bx^. + C 4 0 



(6-14) 



where 



A - S 3 - Y 2 tJ /2 R2/(n-3) 



/s /S ^ 



2 2 
t / nR^r 1 
a / 2 o 1 



B - 2(x 2 k -S 1 -6 2 k )3 3 - ^73 
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The quadratic equation (6*14) can be solved for x^ q using 
the formula. 




-B± 



/b 2 -4AC 

2A 



if A ^ 0 . 



Since 


x k 


x 2k+l 


can 



D - E 



= ^^ x 2k-l +x 2k+l^ ’ the P rediction interval for 
be written as 

^ ^ D "I - E (6 * 16) 



where 

D = ~ x 2k-l ' B/A 

E = /(B Z -4AC)/A (6-17) 

2 

Note: If A=0, then no interval exists; also when B < 4AC 

E becomes imaginary and the interval is taken to be (-oo ; +oo) . 

In section A, a numerical example is given to verify the 
model developed here, with the help of simulated samples. 



B. NUMERICAL EXAMPLE 

In order to compute the prediction intervals for a future 
value for samples satisfying a linear trend model an APL 
program LINTREND has been developed. The listing of this 
program is given on page 67. Error terms having the correla- 
tion structure (6-3) are given as input and for given values 
of a and g, the program develops the vectors and X£. 



For a sample of size 18, the output of the above program is 
given below. 



a = 1. 0 , 6=0.5 



P = 0.2 , N = 18 , n = 8. 




(2.094, 3.37, 



4.056, 3.7595, 5.984, 7.848, 
7.821, 8 , 388) T 



x 2V = (2.94, 3.983, 
8x1 



4.913, 2.685, 5.376, 7.849, 9.017, 
8 . 416) T 



1 



1 



1 2 

1 3 



The design matrix 




1 



4 

5 

6 



1 7 

1 8 



2.094 

3.37 

4.056 

3.759 

5.984 

7.848 

7.821 

8.388 



For the model (6-4), 



£=(-0.147, -0.597, 1 . 566) T 
3x1 



The coefficients of the quadratic equation (6-14) 

A(X^) 2 + B(X£) + C = 0 , 
are given by 
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A = 1.443 
B = 35.801 
C = -72.382 

and the values of D and E are 
D = 15.955 
E = 28.561 

from the equation (6-17) , 
the lower limit = - 12.606 

the upper limit = 44.516 

The true value of (.10.515) is contained in the 

interval . 
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I. VARIATION OF TYPES GF INTERVALS WITH ROW 



SAMPLE SIZE = 14 

SIGMA = 1.0 

LEVEL OF SIGNIFICANCE = 0.05 



ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PRCB 1 


PR03 2 


CCN REG 


0.100 


0.193 


0 .066 


0.741 


0. 9i7 


C. 576 


0. 95 6 


J.2GC 


C.342 


0.080 


0. 578 


0. 91 5 


C. 775 


0. 953 


0. 50 C 


C. 269 


0. C98 


C.633 


0. 883 


0.736 


0 .949 


0. 70 0 


0. 369 


0.119 


0.512 


0.902 


0 .374 


0.949 


0.90 0 


0.539 


0 .128 


0.333 


0.939 


0. 930 


0. 95 E 



Ii. 


VARIATION OF TYPES 


OF INTERVALS WITH ROW 


SAMPLE 


SIZE = l 2 




SIGMA 


= 1.0 




level 


CF SIGNIFICANCE = 0. 


05 



ROW 


'YPE 1 


TYPE 2 


”Y 0 E 3 


D R0 8 1 


PROB 2 


CON REG 


0 .100 


C.4S2 


0. 048 


C. 46 C 


0. S3 5 


C.438 


0. 941 


C. 30 C 


C. 516 


0. C77 


C.407 


0. 924 


0.714 


0 .939 


0. 500 


0. 483 


0.110 


0.407 


0.909 


0 .818 


0.936 


0.70 0 


0.791 


0 .071 


0.138 


0. 934 


C. 873 


0. 939 


0.90C 


C. 902 


0.C5C 


0. 048 


0. 9<t7 


C. 920 


0. 94 8 



CON REG = CONFIDENCE REGION 
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III. VARIATION OF T Y P ES CF INTERVALS WITH ROW 



SAMPLE SIZE = 30 

SIGMA = 1.0 

LEVEL OF SIGNIFICANCE = 0.05 



ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PROB 1 


FROB 2 


CCN REG 


0.100 


0 .55 A 


0.036 


0.410 


0.960 


G. 722 


C. 96 8 


0.30C 


C. 861 


0. 033 


C. 106 


0. 940 


C. 758 


C. 940 


0. 500 


C. 874 


C. C44 


0.082 


0. 952 


0.341 


0.951 


0.700 


0.914 


0.028 


0.058 


0.951 


0.9 29 


0 .953 


0.900 


0.979 


0 .003 


0.013 


0. 952 


1.000 


0. 953 


IV. 


VARIATION OF T 


Y P ES OF 


INTERVALS 


WITH 


ROW 


SAMPLE 


SIZE 


= 34 










S I C-MA 




= 1. 0 










LEVEL 


OF SIS 


NI FI CANCE 


= 0.05 









ROW 


TYPE 1 


TYPE 2 


T YPE 3 


PROS 1 


PROB 2 


CON PEG 


0.100 


C.5S5 


0. 041 


0. 374 


0. 958 


C. 61 0 


C. 96 0 


0.30C 


C. 68 3 


0. C66 


0.251 


0. 943 


0.773 


0.946 


0.500 


0.718 


0.077 


0.205 


0.937 


0.766 


0.837 


0.700 


0.950 


0 .013 


0.037 


0. 942 


0. 846 


0. 943 


0.90C 


C.992 


0.004 


0. 004 


0. 94 8 


1. 000 


C. 94 8 
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V. VARIATION] 0= t Y PES OF INTERVALS WI T H ROW 



SAMPLE SIZE =38 

SIGMA =1.0 

LEVEL OF SIGNIFICANCE = 0.05 



ROW 


T YPE 1 


T Y PE 2 


type 3 


D R0 B 1 


PR0B 2 


CON REG 


0.100 


C. 703 


0. 030 


0. 267 


0. 949 


C.467 


0. 54 8 


0.30C 


C. 809 


0. C36 


0.155 


0. 954 


0.806 


0 .956 


0.500 


0. 847 


O 

• 

o 
»— * 


0.112 


0.956 


0 .8 54 


0.957 


0.700 


0.96 6 


0 .015 


0 .019 


0.959 


0.933 


C. 955 


0.900 


C.994 


0. C04 


0. 002 


0. 957 


1. 000 


C. 557 


VI 


:. VARIATION OF TYPES OF 


INTERVALS WITH 


ROW 


SAMPLE SIZE 


= +2 










SI GMA 


= 1.0 










LEVEL CF SIGNIFICANCE 


II 

O 

• 

o 

U1 








ROW 


TYPE 1 


TYPE 2 


type 3 


PR0B 1 


PR0B 2 


CON REG 


G. IOC 


C. 573 


0. C47 


C.38C 


0.932 


0.723 


0.948 


0. 300 


0.719 


0. 044 


0.237 


0.947 


0.773 


0.952 


0.500 


0.853 


0 .047 


0.10 0 


0.943 


0. 894 


0. 546 


0.70C 


C.980 


0.C07 


0.0 13 


0. S48 


C. 857 


C. 54 8 


C. 90C 


C. 998 


0. C01 


0. 001 


0. 947 


1.000 


0 .547 
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VII. VARIATION OF TYPES OF INTERVALS U I TH ROW 



SAMPLE SIZE = 14- 

SIGMA = 1.0 

LEVEL CF SIGNIFICANCE = 0.10 



ROW 


TYPE I 


TYPE 2 


TY° E 3 


PROB 1 


PROB 2 


CON REG 


C.10C 


C. 269 


0. Ill 


C. 62 C 


0. 87 0 


0.468 


0.906 


0.300 


0.441 


0. 108 


0.451 


0.846 


0 .667 


0.89 6 


0.500 


C.380 


0 . 138 


0.482 


0.858 


C. 732 


C. 9 C 9 


0.700 


0.501 


0. 138 


C.361 


0. 85 8 


C. 790 


C. 900 


C.90C 


C. 666 


0. 114 


C. 22 C 


0. 883 


0.868 


0 .907 


VII I 


.VARIATION 0= 


TYPES CF 


I NTERVALS 


i WITH 


RCW 



SAMPLE SIZE =22 

SI G VA = 1 .0 

LEVEL OF SIGNIFICANCE = 0.10 



ROW 


TYPE 1 


TYPE 2 


T Y° E 3 


PRG3 1 


PRCB 2 


CON REG 


0. 100 


0. 524 


0.077 


0.399 


0.876 


0 .468 


0.894 


0.300 


0.599 


0 .101 


0.300 


0.373 


0. 564 


0. 880 


0.500 


0.617 


0.122 


0.261 


0. 861 


C.754 


0. 884 


0.70C 


C. 853 


0. C66 


0. 081 


0. 88 5 


0.773 


0 .837 


0.900 


0.941 


0. 023 


0.031 


0.991 


0 .786 


0 .89 1 
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IX. VARIATION 0 : TYPES CF INTtRVAlS WITH ROW 



SAMPLE SIZE = 30 

SI GMA =1.0 

LEVEL OF SIGNIFICANCE = 0.10 

















ROW 


TYPE 1 


TY®E 2 


TYPE 3 


PROB 1 


FRQB 2 


CON REG 


0. 100 


0.595 


C.070 


0.33 5 


0.899 


0.571 


0.910 


0.300 


0.906 


0 .025 


0.069 


0.89 1 


0. 560 


C. 65 C 


0.500 


C.935 


0.C24 


0. 043 


0. 902 


C. 875 


C.506 


C. 70C 


C. 951 


0. C22 


0.027 


0.902 


0.909 


0 .905 


0.900 


0.989 


0. 005 


0.005 


0.89 7 


1 .000 


0.898 


X . 


VARIATION OF 


T YRES OF 


INT ERV ALS 


W ITH 


ROW 


SAMPLE 


SIZE 


= 1 4 











SIGMA = 1 .0 

LEVEL OF SIGNIFICANCE = 0.02 



ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PR CB 1 


PROB 2 


CON REG 


0.100 


0.160 


0 .028 


0.812 


0. 969 


0. 571 


0. 583 


0.300 


0.254 


0. C42 


0.704 


0. 949 


C. 738 


C. 576 


C. 50 C 


C. 169 


0. C63 


C. 76 8 


0.929 


0.825 


0 .577 


3.70 0 


0.209 


0.088 


0.703 


0.933 


0 .932 


0.580 


0 .900 


0.3 81 


0 . 121 


0.498 


0.961 


C.975 


C. 58 2 
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XI. VARIATION CF 'YCES OF INTERVALS WITH ROW 



SAMPLE SIZE = 22 

SIC-MA = I . 0 

LEVEL OF SIGNIFICANCE = 0.02 



ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PRGB 1 


FROB 2 


CCN REG 


0.100 


0.464 


0 .026 


0.510 


0. 970 


0.731 


0. 979 


0.30 0 


C.427 


0. 054 


0. 519 


0. 96 7 


C. 870 


G. 979 


G.50C 


C. 338 


0. C76 


0.586 


0. 938 


0.987 


0.978 


0.700 


0.691 


0. 07} 


0.236 


0 .965 


0 .973 


0.974 


0.900 


0.82 5 


0.073 


0.10 2 


0.979 


0. 986 


0. 982 



XII. 

sample 

S I 

LE VEL 


VARIATION OF TY 0 ES OF 

SIZE =30 

= 1.0 

OF SIGNIFICANCE = 0.02 


INTERVALS 


WITH 


ROW 


ROW 


T YPE 1 


"YFE 2 


TYPE 3 


PROS 1 


PRO 8 2 


CON REG 


0.100 


0.525 


0. C19 


0.456 


0. 983 


C. 737 


C. 986 


0.30C 


C. 830 


0. C22 


0. 148 


0. 972 


0.773 


0 .972 


0.500 


0.810 


0. 032 


0.158 


0.978 


0.93 8 


0.98 0 


0.70 0 


0.846 


0 .051 


0.10 3 


0.979 


C. 961 


C. 9 8 C 


0 .90 C 


C.961 


0. 015 


0. 0 24 


0. 9 77 


1. 000 


0. 97 8 
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X III .V AR IAT ION n- -YPES OF INTERVALS WITH ROW 



SAMPLE SIZE = l 2 

S IG^A =3.0 

LEVEL OF SIGNIFICANCE = 0.05 



TYPE I type 2 TY°E 3 PROS 1 PROS 2 CON REG 



ROW 



0.100 
C. 30 C 
0.500 
0 .700 
0 .90 0 



0.067 
C. 179 
0.44 5 
0.736 
C. 901 



0. C78 
0. 125 
0.117 
0 .083 
0. 051 



C. 855 
C. 696 
0.438 
0.176 
0. 048 



0. 851 
0. 855 
0.899 
0.938 
0. 948 



C. 500 
G .704 
0 .829 
C.398 
C. 941 



0.951 
0 .93 7 
0.935 
C. 945 
C.950 



XIV. VARIATION OF ’YPES OF INTERVALS WITH ROW 



SAMPLE S IZE = 26 

SIGMA =3.0 

LEVEL CF SIGNIFICANCE = 0.05 



ROW TY° E 1 TYPE 2 TYPE 3 PROS 1 PROS 2 CON REG 



C. IOC 
0. 300 
0.50 0 
0.70C 
0. 90 C 



C. 652 
0.255 
0.561 
C. 841 
C. 962 



0. C43 
0.106 
0 .110 
0. 062 
0. C 16 



0.305 
0.639 
0.329 
0. 09 7 
0 . 022 



0. 936 
0.929 
0.933 
0. 941 
0. 946 



0.651 
0 .693 
C.855 
C. 903 
0.938 



0 .943 
C.950 
C. 949 
0. 944 
0.947 
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XV. VARIATION 0 : 1 TYPES OF INTERVALS R I TH ROW 



SAMPLE SIZE =33 

SI GMA =3.0 

LEVEL CF SIGNIFICANCE = 0.05 



ROW 



C. IOC 
0.300 
0.50 0 
0.70C 
0.90C 



TYPE 1 



C. 075 
0.374 
0.674 
C. 897 
C. 979 



TYPE 2 



C. C71 
0. Ill 
0 .091 
0. 034 
0. C07 



T Y° E 3 



C. 854 
0.515 
0.235 
C.069 
C. 01 4 



PR OB 1 



0. 880 
3.922 
0.942 
0. 95 C 
0. 952 



PROB 2 



0.620 
0 .721 
C. 857 
C. 912 
1 .000 



CON REG 



0 .964 
0.940 
C. 948 
0. 952 
0.953 



XVI . VARI ATI ON OF TYPES CF INTERVALS WI T H °0‘W 

SAMPLE SIZE = 22 

SIGVA = 5.0 

LEVEL CF SIGNIFICANCE = 0.05 



ROW 


TY° E 1 


TY D E 2 


TY°E 3 


PROB 1 


PROB 2 


COM REG 


0. 100 


0.067 


0.078 


0.855 


0.851 


0 .500 


0.951 


0.300 


0.182 


0 .127 


0.69 1 


0.863 


C. 717 


C. 939 


0.50 0 


C.444 


0. 115 


0. 441 


0. 901 


C. 317 


0. 93 5 


0. 70 C 


C. 716 


0. 102 


0. 182 


0. 936 


0.892 


0.943 


0. 900 


0.901 


0. 052 


0.047 


0.948 


0 .942 


0.950 
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-XVII. VAR I A TI ON 0= TYPES GF INTERVALS WI T H ROW 



SAMPLE SIZE = 26 

SIGMA =5.0 

LEVEL OF SIGNIFICANCE = 0.05 



ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PRQB 1 


FROB 2 


CGN REG 


0. 100 


0.193 


0.097 


0.710 


0.896 


0.619 


0.543 


0 .30 0 


0.255 


0.107 


0.638 


0.929 


C. 701 


C. 550 


G.50C 


C.56 0 


0. Ill 


0. 329 


0. 938 


C. 856 


0. 549 


0. 70 C 


C. 841 


0. C 62 


0. 097 


0. 941 


0.903 


0 .944 


0. 900 


0.962 


0.016 


0.022 


0.94 7 


0.938 


0.548 



XVI I I. VARIA T IQN OF TYPES CF 

SAMPLE SIZE = 30 

SIGMA = 5 .0 

LEVEL OF SIGNIFICANCE = 0.05 


INTERVALS WITH ROW 


ROW 


TYPE 1 


TYPE 2 


TYPE 3 


PRCB 1 


FF03 2 


CCN REG 


0.100 


0.075 


0.075 


0 .850 


0.893 


C. 640 


C. 565 


0 .300 


G .35 0 


0. 110 


0. 540 


0. 920 


C. 7 09 


C. 54 0 


0. 50 C 


C. 664 


0. C97 


C.235 


0. 950 


0.845 


0.952 


0.700 


0.897 


0. 034 


0.069 


0.950 


0 .912 


0.952 


0.900 


0.979 


0 .007 


0.014 


0. 952 


1.000 


0. 553 
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XIX. V ARI A~ I ON GF ""YPES CF INTERVALS WITH N 



CORR.COEFF. = 0.300 

SIGMA = 1.0 

LEVEL GF SIGNIFICANCE = 0.05 



N 


TYPE 1 


TYPE 2 


TYPE 3 


PRCB 1 


FROB 2 


CCN REG 


4 


C .3 00 


0 .040 


0.660 


0. 867 


C.600 


0. 944 


6 


C.342 


0. C80 


0. 578 


0. 91 5 


C. 775 


0. 953 


8 


C. 58 7 


0. C57 


C.356 


0. 942 


0.737 


0.951 


10 


0.516 


0.077 


0.407 


0.92<t 


0.714 


0.939 


1 2 


0.259 


0 .106 


0 .635 


0. 946 


C. 698 


0. 554 


14 


C.861 


0. C33 


C. 106 


0. 940 


C. 758 


C. 54 C 


XX 


. VAR IATION OF TYPES OF 


I NTERVAL S 


WITH 


N 


CGPR. 


COEFF. 


= 0 .500 








SIGMA 




= 1.0 










LEVEL 


CF SIGNIFICANCE 


= 0.05 









N 


T Y 0 E 1 


TYRE 2 


TYPE 3 


PRGB 1 


PRGB 2 


CON REG 


4 


C. 232 


0. C49 


0.719 


0.353 


0.694 


0.951 


6 


0. 269 


0.098 


0.633 


0.388 


0 .786 


0.949 


8 


0.540 


0 .083 


0.377 


0. 935 


C.319 


C. 9 5 C 


10 


C . *83 


0. 110 


C.407 


0. 909 


C. 818 


C. 536 


12 


C. 561 


0. 109 


C.33C 


0. 941 


0.853 


0.95 1 


14 


0.874 


0. C44 


0.082 


0.952 


0 .841 


0.95 1 
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XXI . VARI A TI ON 0 : TYPES CF INTERVALS WI T H N 



CORR .COEFF • = 0.700 

SI GMA =1.0 

LEVEL CF SIGNIFICANCE = 0.05 



N 


TYPE 1 


TYPE 2 


TYPE 3 


PROB 1 


PROS 2 


CON REG 


4 


0.308 


0. 07> 


0.620 


0 .89o 


0 .861 


0.958 


6 


0.36 9 


0 .119 


0.512 


0.90 2 


C. 874 


C. 949 


a 


C.5 86 


0. 095 


C. 319 


0. 920 


C. 884 


C. 942 


10 


C. 791 


0. C71 


0. 138 


0. 934 


C.373 


0.939 


12 


0.840 


0.065 


0.095 


0.940 


0 .908 


0.944 


14 


C. 914 


0 .028 


0.058 


0.951 


C.9 29 


C. 953 


XXII .VARIATION OF T YP£S OF 

CORR.CCEFF. = 0.300 

SIGMA = 3.0 

LEVEL CF SIGNIFICANCE = 0.05 


INTERVALS 


WITH 


M 



N “YPc L T YPE 2 TY D E 3 PROS 1 PRO B 2 CON REG 



4 


C. 067 


0. 050 


0. 883 


0. 642 


C. 580 


0. 955 


6 


C. 090 


0. C93 


0.817 


0.689 


0.710 


0.945 


8 


0.262 


0. 085 


0.653 


0.924 


0.718 


0.956 


10 


0.179 


0 .125 


0.696 


0. 855 


C. 704 


C. 937 


12 


C.255 


0. 106 


C.639 


0. 929 


C. 698 


C. 950 


14 


C. 374 


0. Ill 


0. 5 1 5 


0.922 


0.721 


0.94 0 
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XX III. VAR 1 4 TI ON OF TYPES OF INTERVALS WITH N 



CORR.CCEFF. = 0.500 

SI G^A =3.0 

LEVEL CF SIGNIFICANCE = 0.05 



N 


TYPE 1 


TYPE 2 


T Y° E 3 


PR OB 1 


PR08 2 


CON REG 


4 


C. 108 


0. C 6 3 


C.824 


0. 713 


0.750 


0.952 


6 


0.226 


0. 105 


0.669 


0.863 


0 .781 


0.946 


3 


0.376 


0. 113 


0.506 


0.923 


C. 822 


C. 95 C 


10 


C.445 


0. 117 


0. 438 


0. 899 


C. 829 


C. 935 


12 


C. 561 


0. 110 


C.329 


0. 93 8 


0.855 


0.949 


14 


0.674 


0. 091 


0.2 3 5 


0.942 


0.857 


0.948 



XXIV. VARIATION OF TYPES OF 

CORR.COEFF. = 0. 70C 

SI G P* A =3.0 

LEVEL OF SIGNIFICANCE = 0.05 


INTERVALS 


WITH 


N 


N 


TYPE 1 


TYPE 2 


TYPE 3 


PRGB 1 


PROB 2 


CCN REG 


4 


0.216 


0 .077 


0.707 


0.829 


C. 870 


C. 953 


6 


C .354 


0. 122 


0. 524 


0. 898 


C. 835 


C. 95 0 


8 


C. 557 


0. 118 


C.325 


0. 92 6 


0.864 


0.94 3 


10 


0. 736 


0.088 


0.176 


0.933 


0 .893 


0.945 


12 


0.841 


0 .0 62 


0.097 


0. 941 


C.903 


0. 944 


14 


0.397 


0. C34 


C. 069 


0. 9 50 


C. 912 


C. 952 
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XXV. VARIATION OF TYPES OF INTERVALS WITH N 



CORR.CCEFF . = 3.300 

S I G V A = 5.0 

LEVEL OF SIGNIFICANCE = 0.05 

















N 


TYPE 1 


T YPE 2 


TY°E 3 


PROS 1 


PROB 2 


CON REG 


4 


C.06 8 


0.C51 


0.881 


0. 632 


C. 588 


0. 954 


6 


C. 088 


0. C94 


C. 818 


0.682 


0.713 


0.945 


8 


0.138 


0. 093 


0.769 


0.884 


0 .710 


0.957 


10 


0.182 


0 .127 


0.69 1 


0 . d6 3 


C. 717 


C. 939 


12 


C .25 5 


0. 107 


0. 628 


0. 929 


C. 701 


C. 950 


14 


C. 350 


0. 110 


C.540 


0. 92 0 


0.709 


0.940 


XXVI 


. VARIATION OF TYPES CF 


I NTERVAL 


s with 


N 


CORR . 


COEFF. 


= 0. 50C 








SIGMA 




= 5.0 










LEV EL 


CF S IGNIF ICANCE 


= 0.05 









N 


TY D E 1 


TYPE 2 


TYPE 3 


PROB 1 


PROB 2 


CON REG 


4 


0.107 


0. 069 


0.824 


0 .701 


0 .739 


0.950 


6 


0.215 


0 .109 


0 .676 


0.860 


C. 761 


C. 944 


8 


C .363 


0. 126 


0. 511 


0. 92 8 


C. 81 0 


0. 95 0 


10 


C. 444 


0. 115 


C.441 


0. 901 


0.817 


0.935 


12 


0. 560 


0. Ill 


0.329 


0.938 


0.856 


0.949 


14 


0.664 


0 .097 


0.239 


0. 950 


0.845 


C. 952 
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XX VI I .VAR I ATI CN CF T Y FES OF INTERVALS WITH N 



CORR. COEFF . = 0. 900 

SI G^A = 5 .0 

LEVEL OF SIGNIFICANCE = 0.05 



N 


TYPE 1 


TYPE 2 


TYPE 3 


PRQB 1 


FROB 2 


CON REG 


4 


0.237 


0 .103 


0.66 C 


0.873 


C.922 


0. 562 


6 


0.537 


0. 129 


C. 334 


0. 93 7 


C. 930 


C. 557 


8 


C. 777 


C. C80 


0. 143 


0. 943 


0.888 


0 .947 


10 


0.901 


0. 05? 


0.047 


0.948 


0 .942 


0.550 


12 


C .962 


0 .016 


0.022 


0.947 


C. 938 


C. 546 


14 


C.979 


0. 007 


0. 014 


0. 952 


1. 000 


0. 553 


XXVIII 


.VARIATION OF 1 


NTERVALS 


WITH S I G W A 




S AM°le 


S IZE 


= 14 










CORR.COEFF. 


= 0.700 








LEVEL 


CF SIGNIFICANCE 


= 0 .05 








s igma 


TYPE 1 


TYPE 2 


TYPE 3 


PROB 1 


PROB 2 


CON PEG 


1 


C. 369 


0. 119 


C.512 


0. 902 


0.374 


0.549 


2 


0.352 


0. 123 


0.525 


0.901 


0 .886 


0.951 


3 


0.354 


0.122 


0.524 


0.898 


C. 885 


C. 550 


4 


C .354 


0. 121 


0. 525 


0. 898 


C. 884 


C. 55 0 


5 


C. 353 


0. 122 


0.52 5 


0. 398 


0.885 


0.950 
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XXIX. VARIATION OF INTERVALS WITH SIC-MA 



SAMPLE SIZE =18 

CQRR.CCEFF. = 0.500 

LEVEL OF SIGNIFICANCE = 0.05 



SIGMA 


TYPE 1 


TYPE 2 


TYPE 3 


PROS 1 


FRO 9 2 


CCN REG 


1 


C.540 


0.083 


0.37 7 


0.935 


C.819 


C. 55 C 


2 


C.3 56 


0. 116 


0. 488 


0. 919 


C. 819 


0. 54 7 


3 


C. 376 


0. 118 


C. 506 


0.923 


0.322 


0.950 


4 


0.372 


0. 118 


0.51 0 


0.922 


0.822 


0.550 


5 


0.363 


0 .126 


0.511 


0. 928 


C.8 10 


C. 55C 



XXX 


.VARIATION OF 


IN TE D V AL3 


WITH SIGMA 




SAMPLE 


SIZE 


=2Z 










CORO .C9EFF. 


= 0. 


500 








LEVEL 


OF SI 3 


NIFICANCE = 0.05 








SIGMA 


T YPE 1 


TYPE 2 


T Y°E 3 


OR08 1 


°R0 9 2 


CON REG 


1 


C.483 


0. 110 


0. 40 7 


0. 505 


C. 818 


C. 53 6 


2 


C. 446 


0. 115 


C.435 


0. 901 


0.826 


0.936 


3 


0.445 


0. 117 


0.43 8 


0.899 


0.828 


0.535 


4 


0.446 


0.115 


0.439 


0. 899 


0.326 


0. 535 


5 


0.444 


0. 115 


0.44 1 


0. 901 


C. 817 


C. 535 
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PROB. OF TYPES OF INT . 



1. TYPES OF INTERVALS VS p 




SAMPLE SIZE = 30 

SIGMA = 1.0 

SIGNIFICANCE LEVEL = 0.05 

1, 2 and 3 on the curves indicate the types of Intervals. 
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PROB. OF TYPES OF INT. 



2 . 



TYPES OF INTERVALS VS N 




N - 

CORR. COEFF = 0.5 

SIGMA = 3.0 

SIGNIFICANCE LEVEL = 0.05 
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PROB. OF TYPES OF INT . 



3. TYPES OF INTERVALS VS N 




CORR. COEFF. = 0,7 

SIGMA = 1.0 

SIGNIFICANCE LEVEL = 0,05 
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PROB. OF TYPES OF 1NT . 



4. TYPES OF INTERVALS VS p 




CORR. COEFF. p - 

SAMPLE SIZE = 8 

SIGMA = 3.0 

SIGNIFICANCE LEVEL = 0.05 
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) 



PROB. OF TYPE 1 INT. 



5. PROB. OF TYPE 1 INT. VS N 
(VARIATION WITH a) 




CORR. COEFF. = 0.9 

SIGMA = 1.0 
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C PROGRAM TO CALCULATE PROBABILITY CF TYPE 1, T YPE 2 »TY 3 E 

C 3, T YPE A INTERVALS F°R SIMULATED SAMPLES. 

C PROGRAMMER T S MURTHY SEP 1979. 

Cl MENS ION Z(5 5) ,X ( 55 ),V( 100,5 5),XI < 50) ,YI( 50 ) ,S< 55) , 
1IC( 5), IP ( 5), I»C( 5) , STA T( 10, 7) , VS ( 1 0 , 10 , 7) , I VS ( 1 C) 

CALL OVF LOW 
INDEX =1 
SIGMA=1.0 

1 READ( 5,2) K,T 
WRITE ( 6 ,2 ) K » T 

2 FDR maT( IX, 12, 2X =F 5. 3) 

IF { K .EQ. 0 ) GO TO 460 
ROW =0.10 
DO 30G 13=1, 5 
DO 10 J=l,5 
3 = ( 1 - ROW* *2 ) **. 5 

C SIMULATION OF SAMPLES 

I SEED= 12345 
10 I PC ( J ) =0 

CD 250 I 4=1, 10 

DO 5C M=l, 100 

CALL SNORM{ISEED,Z,K) 

X ( 1 ) =SIGMA*Z(1) 

CD 30 J=2 ,K 

30 X ( J )=ROW*X( J- l)f 3*S IGMA*Z( J) 

DO 40 L=1 , K 
40 V(M,L)=X( L) 

50 CONTINUE 

DO 60 11=1,5 
6 0 I R ( 1 1 ) =0 

CD 200 1=1, IOC 
DO 70 J = 1 , K 
70 S(J)=V(I,J) 

K= K -5 
KK=K/2 
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80 



DO EC L=1 i KK 
Y I ( L) «S ( 2 * L J 
N=KK-1 

DO c 0 LL=1,N 
90 >1 ( LL)=< S(2*LL-1) + S (2*LL+1)) /2 .0 

XSUN=0. 0 
YSU^=0 .0 
SXX=C. 0 
SX Y =0 .0 
$YY=C. 0 
DO ICO KL= 1 i 5 
100 IC(KL)=0 

DO 110 M= 1» N 
YSUN*=YSUM+YI ( V) 

XSU*=XSU^ +X I ( M J 
110 CONTINUE 

X3 = > SUM /N 
Y9=YSUM/N 
DO 120 M= It N 
S X X= S XX+ (XI ( M ) -X 6 ) * * 2 
SYY = SYY+(YI(M)-Y3)**2 
SXY=SXY+( XI (V )-X3 ) * ( YI ( M) -Y3 ) 

120 CONTINUE 

V3ES=(SYY-( (S XY**2 )/SX X ) )/ (N-2 ) 

EH= SXY/SXX 

AH=YB-BH*X3 

SC= (7£~2 ) * V R ES 

A=(EH**2 )-SS/S<X 

°= S ( 2 *KK )- YB 

3=— ( 2*3 H* P ) 

C= (°**2 J- (S S* (N+l ) )/N 

F=(eh**2 )*SXX/SS 

SSQ = 1.0-N*( P* *2) / ( SS* ( N+l ) ) 

IF (F . L T . SS C) GO TO 150 
IF ( F .EQ.l .0 ) GO TO 500 
D= ( 2 «*XB ) -S( K-l) +2. *BH*P/A 
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E=( 2. / A)* (( * S$*{ P**2) /SXX)+(N+1)*A*$S/N)**.5) 

P I L =D-E 
PIR=0+E 
P VA L = S( K+ 1) 

IF ( F . G T * 1.0) G TO 130 
1C (2 )= 1C ( 2 ) +1 

1 F { PVAL.LE.PIR .OR. PVAL.GE.PIL) IC( 5) =IC (5 ) + l 
GO ~Q 160 

130 IC(1 )=IC(1 )+l 

IF(DIL.LE.PVAL.AMD.PIR.GE.PVAL) IC(4)=IC(4)+1 
GO TO 160 

150 JC(3 ) = IC ( 3 ) +1 
160 00 17 G J= 1, 5 

170 IR ( J)=IR( J)+I C(J) 

K = K + 5 

200 CONTINUE 

CO 220 J = l» 5 

22 C I PC < J)=IPC( J)+IR( J) 

250 CONTINUE 

C PRINT STATISTICS 

STAT( IB , 1 ) = P0 Vn 
STAT( IB » 2) =1 PCI1) /I CCO. C 
ST A” ( I B t3 } =1 PC (2 ) /I 000.0 
ST AT ( 1 8, 4 ) = IP C(3 )/1000.0 
I F ( STAT( IB ,2) .EQ.O. C) GO TO 275 
STATI IB ,5)=I PC (4) /( ST AT (IB,2)*1000.0) 

GO TC 280 

275 ST AT ( IB» 5 )= I D C( 4 ) 

280 IF ( STAT( 13,3) .EQ. 0.0) GO TO 285 

STAT(IB t 6)=IPC(5) /( ST AT (IB»3)*1000.C) 

GG TC 290 

285 ST AT ( IB, 6 )= I»C<5 ) 

290 SIG^STATI IB ,2 )*ST4T( IB , 5) + STAT(I B,3)*STAT( IB,6) 
1 +STATI IB ,4) 

ST A T ( IB, 7 )=S IG 
RO W=R O W +0 .2 
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300 CONTINUE 

DO 305 I Q=1 ,5 
CO 305 JQ=1,7 

VS (INDEX, 10 » JO) = ST4 T( IQ,JQ) 

305 C JN T I NU E 

I VS ( INDEX )=N 
INDE X=I NDE X+- 1 
I S Z = K-5 

WRITE(6, 310) ISZ. SIGMA, N 
UR ITE ( 6,325 ) 

UR I 7E( 6,3 50) ( (STA T ( K,U »L=1 ,7) ,K = 1 ,5) 



310 


FORMAT (IX 


SAMPLE SIZE = ' 


, 15 , • S I GM A = 


'» F5.0, 


« 




IN = • , 15, /) 










325 


FORMAT! IX,' 


10W TYPE 


1 T YPE 2 


TYPE 3 






lFROe.l FRCB 


.2 CON. REG 


' ,/ ,70 ( «- ' ) ) 






350 


FORMAT ( 7 ( F9 


.3, 2X ),/ ) 









WRITE (6, 485) 

GO TO 1 

460 C C R F =0 . 1 

IN C EX= INDEX-1 
CO 470 J = 1 , 5 

WRITE (6,472) CORR, SIGMA 
UR I T E (6 ,475 ) 

00 471 1=1, INDEX 

UR I TE( e ,4 80 ) ( I V S ( I ) , (VS (I , J ,K) ,K = 2 ,7J J 
47 L CONTINUE 

COR F= CORR +0 .2 
470 CONTINUE 

472 FORMAT ( * CORR.COEFF. = ',F5.3,' SIGMA = *,F5.0,//) 

475 FORMAT ( ' SAMPLE SIZE T YOE 1 TY ^ E 2 TY°E 3 ?R0 

16.1 drop . 2 CON. PEG ',/,70('-')) 

480 FORMATtl 5 , 5X , 6( F 8. 3 ,2 X) ,/) 

500 STOF 
ENC 
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APL PROGRAM LINTREND 



V LIMTREMD X 

Cl] x*-l + <x + 0.5x ( > 

C2] 2 

C3] < h -2 ) fx 

C 4 ] Kf(fX)x 0.5 

C5] J < x <- x 0 

C6] m*-1 

C7] -»1 1 X l K < M<-M+l 

C8] X1<-(X|- ( (2xM)-3) 3 + xc ( (2 xm)-1)])x0.5 

C9] xif-xi ,xi 
CIO] -»7 

C 1 1 ] ti<-xc2x \ K3 
C12] M <- K — 1 
C 13] YlfNfYI 

C14] 1 the vector xi is >jxi 

C15] 1 THE vector yi is • ;n 

Cl 6] 1 the value or xp is • ;xr 

C17] CK-Mfi 
C18] C2<-\H 

C 19] C3<-xi 

C20] fMT<-(3,n) f ci,c2,C3 

C21] DMfUDMT 
C22] SfDMT+,xDM 
C23] Si<-gs 

C24] BHf (SI+, Xl'MT) + , X'l'I 

C25] SHSf(+/( ( f I ) * 2 > >-( (*5TI+, X»M) + ,XBH) 

C26] «0<-sid ? 1] + (M + 1 ) X ( (2x51 Cl J2] >+<"+1 ) xsic2»23 ) 

C27] ao<-«0 + 1 

C28] «ir(2xSIC 1 53] )+2xsiC2?3]x (M+l ) 

C29] A 2<- s i C3 i 3] 

C30] CF«-4xSHSr(M-3) 

C31 ] a«-(&hc3]*2)-( a 2xcf) 

C32] b<-(2xbhc3]xW4-(XC2xK]-bhc1]-bhc2]xi<) )-Aixcr 

C33] Cf<w*2)-fi0xCF 

C34] ->im x \ (B*2)14x«xc 

C35] e<-( (B*2)-4xa x c) *o.5 

C36] 

C37] eh-~1x;:c (2xk >-13-B4-a 

C38] ' the left limit is 1 $d-e 

C39] ' the right limit is • ;d+e 

C 40] ->0 

C41] I M J 1 NO INTERVAL EXISTS 1 

V 
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APL PROGRAM PREDICT 



C 1 3 
C23 
C 3 3 
C43 
C53 
C 6 3 
C73 
C83 
C93 
C 1 0 3 
Cl 13 
C 1 23 
C 133 
C 1 4 3 
C 153 
C 163 
C173 
C 183 
C 1 9 3 
C203 
C21 3 
C223 
C233 
C243 
C253 
C26 3 
C273 



p P: e r< i c r x ; a ; b ; c ; t < ; i :;: ; f ; g ; m j p ; s s«;s s x;s 

n X is ASSUMED TO BE A VECTOR OF EVEN 
K <- < f « ) X 0 . 5 

M«-l 

•4 9 X \ K < Mf-M-H 

( (2XMJ-3) n + i.'C ( (2XM)-1 ) 3 ) xO . 5 

xi«-xi f xi 
-45 

T 1 p *■' C 2 x \ k ] 

MfK-1 
YifMfyi 
'T'Bf ( + /YI ) tM 
KBf (+/KI ) pN 

SXt>+/( (XI-XB) X ( YX-YB) ) 

SXXp+/( (XI-XB) *2) 

5YYf+/ ( ( YX-YB ) *2 ) 

BHf SXY^SXX 
AH<-YB-BHxXB 



x f cxy ) syt ; x b ; tb 

NO OF ENTRIES 



VRESf(SYY+(-lx ( ( SXY*2)rSXX) ) ) ~U~2 

e<- 12.706 4.303 3.182 2.776 2.571 2.447 2.365 2.306 

2.201 2.179 2.16 2.145 2.131 2.12 2.11 2.101 2.09 
2.08 2.074 2.069 2.064 2.06 2.056 2.052 2.043 2 
2.041 2.04 2.04 2.04 2.035 2.03 2.03 2.03 
2.025 2.025 2.02 2.02 2.02 2.02 2.015 2.01 2.01 
t«-G[;w-23 

A<- ( BH*2 ) -< ( -rSXX )xSf( ( T*2 ) XVRES ) ) 

Bf"lx ( 2 X BH x p 4- ( x C 2 X K 3 - TB ) ) 



[283 c«.(P*2)-sx((H+i)iH) 

C29 3 Ff-<BH* 2 ) xsxx~s 

C303 1 the f statistic is ■ jf 

C31 3 SSQ<-1 — (N-r (N+l ))X(F-«2)45 

[323 -»T1X\F>1 

C333 -4 t 2x i F=i 

C343 -4T3x\f<ssq 

[353 1 the interval is of tyre 2 ’ 

[363 -40 

C373 tu-ho 

[383 ( 2xxb ) -N[ ( 2 x K )-13+2XBHxP-ffl 

[393 £<-(4fl)x2x(Sx(( ( F * 2 ) *S3<x ) + < (N + l)?N)xfl) ) *0.5 

C403 1 THE left limit is 1 ;b-e 



C413 


• THE RIGHT 


LIMIT IS 


i 


JI' + E 


C423 


• THE LENGTH 


OF* INT IS 


i 


f 1 (2XE) 


C433 

[443 


-»o 

T2 : 1 THE I NT 


IS OF TYPE 


4 ' 




C453 
C 46 3 


•*o 

T3 ♦ * THE INT 


IS OF TYPE 


3 


< 



v 
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